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Abstract 

Some very unusual features of the hadron structure functions, obtained in the gener- 
alized QCD sum rules, like the surprisingly strong difference between longitudinally 
and transversally polarized p mesons structure functions and the strong suppres- 
sion of the gluon sea in longitudinally polarized p mesons are discussed. Also the 
problem of exact zero contribution of gluon condensates to pion and longitudinally 
polarized p meson quark distributions is discussed. 



It is well known, that in QCD DGLAP equations predicts the evolution of these 
distributions with Q 2 but not the initial values at some relatively low Q 2 from which this 
evolution starts . To obtain these initial values (inputs) usually any form of quark (valence 



and sea) and gluon distributions with some free parameters are assumed. Then by using 



Q \ DGLAP equations, quark and gluon distributions are calculated at all Q and x and 
compared with the whole set of the experimental data, and in this way the parameters 
are fixed. Evidently, such an approach is not completely satisfactory from theoretical 
^ ■ point of view, because, we can loose some unexpected physical features of distributions, 
so it would be desirable to determine the initial distribution directly from QCD. This 
talk will be devoted to the problem of the determination of quark and gluon distributions 
in hadrons in model independent way in QCD, and we will see, that some unexpected 
features of quark distributions really appear. 

A method to determinate valence quark distribution in QCD sum rules at intermedi- 
ate x was suggested in [T] and generalized in Due to this generalization it became 
possible to calculate valence quark distributions in pion [2] and transversally and longitu- 
dinally polarized p- meson and also significantly improve the result for valence quark 
distributions in proton 0j. 

Let me at first briefly remind the method. We start from consideration of the 4-point 
correlator corresponding to the forward scattering of hadron and electromagnetic currents. 

n = -if d A xd^yd A ze^ x+i ^ z (0\T{j h (x), j el (y), f l (0), j h (z)} |0> (1) 

Here p\ and p 2 are the initial and final momenta carried by hadronic the current j , q 
and q f = q + Pi — P2 are the initial and final momenta carried by virtual photons (Lorentz 
indices are omitted) and t = (pi — p 2 ) 2 = 0. 
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The idea of the approach (in the improved version) is to consider the imaginary part 
(in s-channel) of this four-point correlator U(pi,p2,q,q'). It is supposed that virtualities 
of the photon q 2 ,q' 2 and hadron currents p\-,p\ are large and negative \q 2 \ = \q' 2 \ ^> 
\P1\1 IpII ^ ^c, 2 1 where R c is the confinement radius. It was shown in [Ij that in this case 
the imaginary part in s-channel [s = (pi + q) 2 } of U(pi,p2~, qi, q') is dominated by a small 
distance contribution at not small x, so OPE is applicable. 

To find structure function one should compare dispersion representation of the forward 
scattering amplitude in terms of physical states with those in OPE and use Borel transfor- 
mation. Though we should consider the case of forward scattering, i.e. p\ = P2, but, and 
it is significant point of the method, we should keep pi not equal to p\ in all intermediate 
stages and only in final result take the forward scattering limit. Only in such way it is 
possible to effectively (exponentially) suppress all terms in sum rules expect those which 
are proportional to structure function and achieve good accuracy. Finally equating the 
physical and QCD representations (for detail see [2]) after double borelization (on p\ and 
p 2 ) we found 

Im U° QCD + Power correction = 2ttF 2 (x, Q 2 )g 2 h e h K (2) 

Here ImX]°Q CD is the bare loop contribution, (where continuum contribution is elimi- 
nated) and gh is defined as (0\jh\h) = g^. In what follows, we will suppose Mf = M| = 
2M 2 . 

Let us start with the case of the pion. To find the pion structure function one should 
choose the imaginary part of 4-point correlator with two axial (hadron) and two electro- 
magnetic currents and consider the invariant amplitude at tensor structure, P^PyPxP^/ u, 
where P = (p± + p2)/2, fi,u are vector current indexes and A, a -hadron current indexes. 
I shall briefly note the main points of calculations, for detail see [2]. In QCD part of sum 
rules we take into account the following terms: 

1. Unit operator contribution (bare loop, Fig.l) and leading order (LO) perturbative 
corrections proportional to ln(Q 2 /fi 2 ) , where fi 2 is the normalization point. In what 
follows, the normalization point will be chosen to be equal to the Borel parameter fi 2 = 
M 2 = M 2 = 2M 2 . 

2. Power corrections - higher order terms of OPE. Among them, first of all we should 
account the dimension-4 correction, proportional to the gluon condensate (0\^-G^ u G™ u \0). 
Corresponding diagrams are shown on Fig. 2. But, surprisingly, it was found that these 
diagrams contribution to the sum rule is exactly cancelled after Borel transformation on 
p\ and p\. I want to note, that diagrams in total may not cancel, but only that part 
of them, which have imaginary part both p\ and p\ (i.e. just those which contribute to 
the pion structure function). So we see there is no gluon condensate contribution to the 
sum rule for pion structure function. The fact of exact cancellation of contributions of 10 
diagrams on fig. 2 is surprising, but of course one can threat this as a just accidence. But 
just the same take place for operators of the next dimension! 

There are a large number of loop diagrams, corresponding to d — 6 corrections. First 
of all, there are diagrams which correspond to interaction only with gluon vacuum field, 
i.e. only with external soft gluon lines (see Fig. 3). They are, obviously, proportional to 
(g z f abc G a G^pG ) . But one should also take into account diagrams of Fig. 4 which are 
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proportional to (0 | D p G pu D T G^ | 0) and (0 | G a pv D p D T G a a(5 | 0) . Using the equation of 
motion it can be expressed in the terms of 

[(0 | gfif, | 0)] 2 and (0 | g 3 G^G b up G c p J abc | 0) 

As it was shown in [2] , if one take the sum of all diagrams, then the terms, proportional 
to gluon condensate (0 | g 3 G a tll ,G h lip G c pp f abc | 0) exactly cancel each other just in the same 
way, as for d = 4 operator we have discussed above. I want to note, that there is about 
a hundred different diagrams, so it is hard to believe that such a cancellation is again an 
accidence. 

So we can fixed first very unexpected feature of pion dynamics - gluon vacuum contri- 
bution to pion valence quark distribution is exact zero for dimensions 4,6. I want to add, 
that, as we will see a little later, just the same effect appear for longitudinally polarized 
p meson, so this is not the specific feature of pion. It is very like that here we meet a 
consequence of some symmetry, though today we cannot understand it. 

Of course except the diagrams on fig. 2a,b there are another diagrams of dimension 6, 
directly proportional proportional to the four-quark operators a 2 = ct s (27r) 4 ({0|^^|0)) 2 . I 
will not discuss them in this talk, one can see papers j2], 0- Finally quark distribution 
function (hereafter we will note it as QDF) has the following form: 
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/ s 3 M 2 

XU n {x) = 7T^-FT X I 1 - x ) 



1 + 



'q s (/i 2 ) ■ln(Ql/^ 2 
3vr 



x N ( r )(i _ f-so/M^ _ 47ra,(/x 2 ) ■ 47ra 2 u(x) 
l Jl ' (27r) 4 • 3 7 • 2 6 • M 6 x^l-xf 



(3) 



where N(x) = (1 — x)(l + 4xln(l — x)) — 2x(l — 2x)lnx and u(x) is the fourth degree 
polynomial in x, This function u n (x) may be used as an initial condition at Q 2 = Ql for 
solution of QCD evolution equations. The result is shown on fig. 5 (thick line). 

In the numerical calculations we choose Qq = 2 GeV 2 (for pion) and a 2 = 0.23GeV 6 , 
found recently from r decays analysis -see jS], [HI; [II- The continuum threshold was 
varied in the interval 0.8 < s < 1.2GeV 2 and it was found, that the results depend only 
slightly on it's variation. The analysis of the sum rule shows, that it is fulfilled in the 
region 0.2 < x < 0.7, where both power correction and continuum contribution are less 
than 30% and the stability in the Borel mass parameter M 2 dependence in the region 
OAGeV 2 < M 2 < 0.7GeV 2 is good. The accuracy of this result was estimated in [2] to be 
20 — 30%. Comparison with various experimental result (see 0) lead to conclusion, that 
agreement with experiment is good. I want to note, that this result for QDF is based 
only on the QCD sum rules. 

It is possible also to estimate the second moments of QDF. Assume, that u^x) ~ 
1/t/x at small x < 0.15 according to the Regge behavior and u n (x) ~ (1 — x) 2 at large 
x > 0.7 according to quark counting rules. Then, matching these functions with our 

result, one may find the numerical values of the first and the second moments of the 

i i 
w-quark distribution: M. 1 = J u n (x)dx ~ 0.84, M.2 = J xu n (x)dx ~ 0.21 where the values 

o o 
various slightly at reasonable changes of QDF behavior at large or small x. The moment 

M.\ has the meaning of the number of u quarks in 7r + and it should be Aii = 1. The 
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deviation of M\ from 1 characterizes the accuracy of our calculation. The moment M.2 has 
the meaning of the pion momentum fraction carried by the valence u quark. Therefore, 
the valence u and d quarks carry about 40% of the total momentum what is close to 
experimental results. 

In the same way one can calculate valence w-quark distribution in the p + meson. The 
choice of hadronic current is evident. j p = U'j^d and the matrix element is (p + | | 0) = 

2 

where m p is the p-meson mass, g p is the p — 7 coupling constant, g 2 /47r = 1.27 and 
e p is the p meson polarization vector. 

In the non-forward amplitude the tensor structure for determination of w-quark distri- 
bution in the longitudinal p meson is to P p P u P a P\, while for transverse p it is —P^PuSxa) 
(seeQ). 

In the case of longitudinal p meson the tensor structure, that is separated is the same 
as in the case of the pion. It was shown in |3] , that w-quark distribution in the longitudinal 
p meson can be found from Eq.(3) by substituting — > m p , f n — > m p /g p . Sum rules for 
u p (x) are satisfied in the wide x region: 0.1 < x < 0.85 with high accuracy (about 10%). 
Fig.5 (curve with squares) presents xu p (x) as a function of x. The values M 2 = 1 GeV 2 
and So = 1.5 GeV 2 , Qq = 4 GeV 2 were chosen. It is clear, that in this case again gluon 
condensate contribution is cancelled. If we construct the moment of QDF for p^ in the 
same way as for pion then we found M\ = 1.06 and M 2 = 0.39 (for one valence quark). 
One should note, that accuracy of this prediction of moment of pi is very high, because 
in this case sum rules for QDF are covered almost all x region, as I have noted before, so 
extrapolation procedure contribution numerically is negligible. 

So we can fix the second very unusual effect: the valence quarks carry about 80% of 
the total momentum of longitudinally polarized p meson, so gluon sea should be strongly 
suppressed (less than 20%). 

Let us now consider the case of transverse p-meson, i.e., the term proportional to 
the structure P p P u 5\ a . The procedure of calculations are the same except one significant 
points. 

1. In contrast to the pion case, the (G a G a ) correction for transversally polarized 
p(Pt), ■ 

2. In contrast to n and p^-meson cases, the terms proportional to (G^G^), as long 
as (0 I g 3 f abc G a pu G b up G c m I 0) are not cancelled for p T . But (0 | g i Gl\ v G b vp G c pp f abc | 0) is 
not well known; so we need to use estimations based there are only some estimates based 
on the instanton model. Result for QDF in transversally polarized p meson is shown on 
Fig.5 for the region 0.2 < x < 0.65, (line with asters). One can find complete analytical 
form and it detail analysis in our paper 3J. The choose of parameters are the same as in 
previous cases. The main sources of inaccuracy are d = 6 gluon condensate (gives about 
20%) and d = 4 gluon condensate (due the inaccuracy of it own value about factor 1.5, 
which lead to uncertainty in QDF about 20%). Accuracy is worse, about 30 — 50% (better 
in the middle of region and worse at the end). From fig.5 one can see, that difference 
between QDF in the longitudinally (p^) and transversally (pr) polarized p mesons is very 
large, many times larger that uncertainties. 

So, and this is third surprise, we can conclude that QDF significantly depend of 
polarization. Such great difference between QDF of pL and pr is very difficult to explain 
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in usual quark model. I want to add, that to check this effect, we calculate the second 
moment of QDF of pi and px in quite another way, considering of 2-point correlator in 
external tensor field. This method was offered in [S] for pion case many years ago (for the 
method itself of the QCD sum rules in the external fields see f5J). In the [TUj we use this 
method for polarized p meson case. The results totally confirm the conclusion that QDF 
of pl and pt are quite differ : the second moment M 2 for one p L is about 84% (compare 
with 78% discussed above) while M 2 for p? is much smaller and close to usual hadron case 
(about 40% - 50%). (More detail this results and comparison of this two methods one 
can found in ^Dj.) That why I think, that QDF indeed strongly depend on polarization, 
at least in p-meson case, and this fact can't be explained in the present quark models of 
the hadrons. 

Summarizing, in this talk I tried to pay your attention to very strange results, which 
appear from pure QCD non-model calculation of quark distribution function and yet have 
no satisfactory physical explanation. These are: 

a) First of all, the total cancellation of diagrams with gluon condensate (for dimensions 
d = 4,6) seems to note to some unknown symmetry. Of course it is interesting to check 
this property also for dimension d = 8, but it is technically very hard and isn't complete 
yet. 

b) Second, such a strong difference between QDF for longitudinally polarized (pi) 
and transversally polarized (pr) , as is shown on fig.5, can not be explained in the usual 
hadron quark models. 

Taking in mind also the third very strange fact that, as we discuss before, for p^ gluon 
sea is strongly suppressed (less than 15%), we can conclude quark dynamic in (p L ) and 
(Pt) should be quite different, in contradiction to usual quark models. From my point 
of view it is extremely interesting to understand this strange features, at least in some 
model way. 

Author is thankful to B.L.Ioffe for useful discussions. 
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Figure 1: diagrams, corresponding to bare loop contribution 




Figure 2: diagrams, corresponding to d = 4 gluon condensate 
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Figure 3: diagrams, corresponding to d = 6 condensate 




Figure 4: diagrams, corresponding to d = 6 condensate 
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Figure 5: Quark distribution function for pion, transversally and longitudinally polarized 
p - solid line, line with asters, and line with squares correspondingly 
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